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Abstract: In this paper, the calorimetric power measurement method for electron cyclotron resonance 
heating system on EAST are presented. This method requires measurement of the water flow through the 
cooling circuits and the input and output water temperatures in each cooling circuit. Usually, the inlet 
water temperature is controlled to be stable to get more accurate results. The influence of the inlet water 
temperature change on the measurement results is analyzed for the first time in this paper. A novel 
temperature calibration method is proposed also. This kind of calibration method is accurate and effective, 


and can be easily implemented. 
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I. INTRODUCTION 


A 140GHz electron cyclotron resonance heating system for EAST (Experimental Advanced 
Superconducting Tokamak) is being built in ASIPP (Institute of Plasma Physics, Chinese Academy of 
Sciences)[1]. This project is designed to inject 4MW/140GHz/100s RF power to EAST. The first two 
gyrotrons have been established already. The two gyrotrons are produced by Gycom and CPI respectively, 
and they have been tested in ASIPP. To get an accurate interpretation of the energy balance of the plasma 
and its response to heating during electron cyclotron resonance heating (ECRH) experiments, a reliable 


estimate of the RF power delivered to the plasma is needed[2]. 


There are many methods to measure the millimeter wave power, such as calorimetric method (which 
may use thermistor sensors, thermoelectric transducers[3], etc.), detection method (which may use diode 
sensors), etc. These methods which use thermistor sensors, thermoelectric sensors and diode sensors can 
only measure small power. They can not be used to measure the output power of the gyrotrons directly 
because the output power is in MW-level. We can use the couplers to couple small amounts of microwave 
energy from the high power to measure the real power. In addition, we can make big dummy loads which 
are cooling by water flow for gyrotrons, then we can measure the output power of the gyrotrons by 
measuring the calories of the dummy load. This method is called the flow calorimetric method, or 


calorimetric method for short in this paper. The calorimetric method is used to measure the output power 
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of gyrotrons in EAST ECRH system. This method is accurate only when the inlet water temperature is 
stable. The influence of the inlet water temperature change on the measurement results is analyzed for 


the first time in this paper. 


In this paper, we begin by discussing the calorimetric power measurements. Next, we analyze the 
calorimetric method when the inlet water temperature is variable. Then, some test results are given. 


Finally, we give the conclusion. 


ll. CALORIMETRIC POWER MEASUREMENTS 


The calorimetric measurements require the measurements of the water flow through each cooling 
circuit and the input and output water temperature in each cooled circuit. We have established a stable 
data acquisition system to acquire the water flows and the temperatures in each circuit[4]. It is based on 
the PXI platform. PXI is a rugged PC-based platform for measurement and automation systems [5]. It is 
widely used in the area of measurement[6, 7], data acquisition[8], control[9, 10], etc. In our system, the 
data is being acquired when the gyrotron is working, being displayed on the monitor in real time, and 


being saved to the data storage server which are ready to be used for successive analysis procedures. 


The architecture of the calorimetric measurement is presented in Fig. 1. In EAST ECRH system, 
temperatures are measured using platinum resistance temperature transmitters to an accuracy of about + 
0.3°C (0°C to 100°C) and a resolution of 0.02°C, water flow is measured using vortex meters to an 
accuracy of about 1% F.S.. These measured data then be send to PXI 6229 data acquisition cards [11] 
through the isolators. PXI 6229 DAQ cards have an accuracy of 3.1mV when the range is -10V to 10 V, 
and having a sensitivity of 97.6uV when the range is -10V to 10V. The isolators have an accuracy of 
10mV. So, when we set the range of DAQ cards to be -10V to 10V, the voltages from OV to 10V are 


proportional to the temperatures from 0°C to 100°C. So, the relative error of temperature measurement 


is, 


H = (2h + ef + ee ) = pan (1) 
T T 
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[e ad li is caused by the DAQ cards. 
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Similarly, the relative error of flow measurement is, 
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where Fy is the preset maximum flow. For instance, Fyu is set as 100 m*/h for Gycom Dummy load. 
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The microwave output power generated by a specific gyrotron can be measured by directing the beam 
to dummy load which is cooled by deionized water. By measuring the inlet water temperatures and the 
outlet water temperatures of dummy load, we can calculate the power absorbed by dummy load [2, 12, 
13]. We can use the same way to get the power absorbed by main window[14], MOU, waveguides, 
waveguide switch, miter bends, and bellows which can be used to estimate the gyrotron output power 
when the rf power is directed to tokamak. The generated power of gyrotron is the sum of those power 
data mentioned above. We also use the same way to get the power absorbed by relief load, and relief 


window. 


Isolators 


Water 
Out 


‘Transmission Line 


MOU 
FIG. 1. Architecture of calorimetric measurements. 
As shown in Fig. 2, tọ is the initial time. The water, whose temperature is Tin(to), run through cross 


section A at time to to point B after time ôt. The temperature of the water is Tout(ty + ôt) at the point 


of B. 
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FIG .2. Schematic of calorimetric measurements. 


We know that the energy rise of the water can be calculated from 
W=C-m-AT (3) 


where C is the specific heat capacity of water, m is the mass of the water, AT is the temperature rise. So, 


from time tọ + ôt to time ty + St + At at point B, the energy rise from point A to point B is, 
Wro+st+at = C =m [Tout(to + ôt + At) = Tin(tg)] (4) 


where m = F - At, F is the mass flow of water in this circuit, it is controlled to be stable in experiments. 


So, F can be seen as a constant here. Similarly, from time to to time ty + At at point B, the energy rise is, 
Wo = C- m: [Tout(ty + At) — Tin (to — &t)] (5) 

So, from time to to time tn, all energy rise from point A to point B is, 

W =$} W; = XC: F At: [Tout(t; + At) — Tin(t; — ôt)] (6) 
where t;4, — t; = At. 
When At > 0, 
W = jim, yr C:F At: [Tout(t; + At) — Tin(t; — 5t)] = jim uC: F- At: [Tout(t;) — 
Tin(t; — 8] =C-F- f,."[Tout(t) — Tin(t — 8t)]dt (7) 
Assume the wave-output duration is T, the average power of the gyrotron over the time span t can be got 
by 


pot a2. > [Tout(t) — Tin(t — ôt)]dt (8) 


T T 


We found that when gyrotron ceases generating rf power, the difference between outlet water temperature 
and inlet water temperature rises still, it comes to the maximal value after a while, then it descends slowly 
to zero within 60 seconds after the pulse is ended, as shown in Fig. 3. This means that energy is absorbed 
and temporarily stored in structural material with possible thermal loss to the outside. The thermal loss 


is assumed to be little enough to be ignored. So, 
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P Tout(t) — Tin(t — ôt)]dt ~ z, Tout(t) — Tin(t — 6t)]dt (9) 
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FIG. 3. The typical signal of temperature difference between output water and input water. The pulse duration is T. 
In the measurement, the water temperature at the inlet is kept as an approximate constant, thus, formula 


(9) can be simplified to, 


_ CF ptn=tot+t+60 


P= 


tn=to+t +601 


j zer; 
[Tout(t) — Tin]dt = a 


Tout(t) — Tin(t)]dt (10) 


T to 
As mentioned before, the absolute error of the measured temperature is about +0.3177°C, so, the 


maximum power error is 


Pyro = E- i +0.6354dt = +. 0.6354: (t + 60) (11) 


T 
But when the pulse is very short, the maximum of Tout(t) — Tin(t) is smaller than 0.6354 °C in some 


circuits. So, the real power, 


Poa = ÊE. Je Tout(t) — Tin(t)]dt < CT. 0.6354: (T + 60) (12) 


T 
That is, the relative error of the power can be up to +100%. In order to solve this problem, we must 
calibrate the temperature difference first. The real temperature difference equals the measured 


temperature difference plus a base temperature, an example is shown in Fig. 4. 
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FIG. 4. The temperature is influenced by calibration. The relative temperature is more accurate than the absolute 
temperature, it is influenced by the resolution of the temperature sensor only. 


In the case of no RF power transmitted to the load, the inlet water temperature should be equal to the 
outlet water temperature. But due to the existence of the measurement error of the temperature sensors, 
the difference between the inlet temperature and the outlet temperature is not zero but an approximate 
fixed value. So we need to calibrate the temperature measurement to eliminate the error. That is, we need 
to find the difference value between the outlet water temperature and the inlet water temperature in the 
absence of RF power. 

There are two methods to get the compensation value. One method is using the digital low pass filter 
to filter out the noise signals of the temperature signals. Then we can get the inlet temperature Tin and 


outlet temperature Tou. So, the compensation value is, 


Te = T; ~ Tout (13) 
In the ideal case, the temperature difference after calibration becomes, 
ATac = Tout — Tin + Te = 0 (14) 
So, ideally, when no RF power is transmitted to dummy load, the measured power will be approximately 
equal to zero. In fact, the ambient temperature is not a constant. Even if the water temperature is 
completely constant, the temperature measured by sensors will always be a little jitter. So, the 
compensation value is not particularly accurate using the method of low pass filter. 

The other method to get the compensation value is using integral formula (10). We give a test pulse 
signal with duration t (taking into account the time measurement error, 10s>t>500ms is appropriate) to 
start the power integral program. Because no RF power transmitted to dummy load, the measured power 
should be zero. But due to the temperature difference between the outlet water and the inlet water, the 


measured power using formula (10) is not zero, but is, 


Paai = fE Toutt) — Tin(t)]dt = Z, AT- (t + 60) (15) 


T to 
Where, AT is the average temperature difference between the inlet and the outlet. An example of the 
waveforms of the calorimetric power measurements in the absence of RF power is shown in Fig. 5. The 
test pulse duration is 0.896s. As we can see, the measured dummy load power is gradually increasing 
with an approximately constant slope. Using the measured power value, the compensation value can be 


calculated by, 


Pcalt 


T, = -AT = -e 
CF(t+60) 


(16) 


All variables in formula (16) are in the International System of Units. If we change the unit of flow signal 


F from ‘kg/s’ to the commonly used unit ‘m*/h’, the compensation value can be calculated by, 


3.6 PcalT 


T, = -AT = 
CF(t+60) 


(17) 


This method is equivalent to find the average temperature difference between the inlet and the outlet. 
This method is more accurate than the first method of low pass filter. And it is easy to implement because 
this calibration method is using power measurement principle to calculate the compensation value. We 
choose this approach to calibrate the temperature measurement. As the ambient temperature and water 
flow are not constant, we need to calibrate the calorimetric power measurement system periodically to 


calculate the latest calibration value to get more accurate measurement results. 


The load power obtained by the 


direct multiplication method. 
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FIG. 5. An example of the waveforms of the calorimetric power measurements of the dummy load. The dummy 


load power was measured in the absence of RF power. 


The temperature difference error is effectively corrected by calibration. But the resolution of 
temperature transmitter can cause measurement error also. Actually, after temperature calibration, the 
power measurement error is caused mostly by the resolution of temperature transmitter and the accuracy 


of the flow meter. Assume, 


tn=to+t +601 


l= Je Tout(t) — Tin(t)|dt (18) 


Then, we have, 


CF 
Preat = — ` Tint (19) 


T 


In the worst case, the relative error of measured power is, 


A O1n(Preal) | |? In(Preal) | + [Erea A z| [ee In(Preal) |= (= "l |Z| + 
Spower ~ + || E ar| + Se ATine : + ac|] = + (F+ F 


“| +) (20) 
c 
where AF , ATint, AT is the absolute error of flow, temperature integral, wave-output duration 


respectively. In our system, as see in formula (2), 


j= pe en 
In our system, the duration measurement error is about 1ms, and less than 2ms in the worst case, 
AT < 2ms = 0.002s (22) 
We set the specific heat capacity of water as, 
C = 4.18 kJ/(kg - °C) (23) 


But actually, the specific heat capacity will change along with the temperature and the pressure, as shown 
in Table 1. And the water temperature is between 0 and 100 °C and the water pressure is between 0.1 


MPa and | MPa. So, the maximum relative error of the specific heat capacity is, 


E __ 4.217-4.18 
cls 4.217 


x 100% = 0.88% (24) 


TABLE 1. Liquid water specific heat capacity at constant pressure along with the temperature and the pressure. 
The unit is kJ/(kg - °C). 


Pressure Temperature (°C) 

ee) g 2 50 700 
0.1 4.217 4.182 4.181 
0.5 4.215 4.181 4.180 4.215 
1.0 4.212 4.179 4.179 4.214 


The typical temperature difference curves are shown in Fig. 6. Assume Tmax is the maximum of the 
difference between final and initial water temperatures in each circuit, t, is the duration when 
temperature difference rise from zero to maximum (It is related to the response time of temperature 
transmitter, the character of load, and the gyrotron pulse duration). As we can see, 

Tint © Max(T, tr ) * Tmax (25) 


Because the temperature resolution error is +0.02 °C, in the worst case, 


AT ine ~ 0.04 (t + 60) (26) 
Hence, 
ATint| _ 0.04:(t+60) 
| Tint max(t,ty )-Tmax (27) 


But due to the randomness of the temperature resolution error, usually, we have, 
AT; np = 0 (28) 


So, the maximum calorimetric measurement error is, 


Spower ~ + (Lt + S280 __ 4° + 0.88) % (29) 


F max(t,tr )'Tmax T 


The probable calorimetric measurement error is, 
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Srower © + (= ++ 0.88) % (30) 
where F is the mass flow of water in each circuit; Tmax is the maximum of the difference between final 
and initial water temperature in each circuit; T is wave-output duration; t, is the duration when 
temperature difference rise from zero to maximum, for instance, t, ~ 3 s for waveguide dummy load of 
CPI gyrotron when the pulse duration is 2 s, t, ~ 14 s for aluminum tank dummy load of CPI gyrotron 
when the pulse duration is 2 s, t, ~ 13 s for dummy load of Gycom gyrotron when the pulse duration is 
5 s. And all parameters are in standard international unit. When the gyrotron pulse duration is bigger than 


1 s, the relative error of calorimetric measurement is less than 10% and less than 5% in most cases. 
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FIG. 6. The temperature difference between outlet and inlet of dummy load cooling circuit. The figure above shows 


the situation when tr is bigger than t. The figure below shows the situation when tr is smaller than t. 


lil. CALORIMETRIC MEASUREMENTS WHEN THE INLET WATER TEMPERATURE IS 


VARIABLE 


When the water temperature at the inlet is variable, the power measurement can be got by formula (9). 
Firstly, the parameter ôt should be measured. Assume the water flow is F, the water volume of the 


dummy load is V, then, we have, 
v 
ôt = z (31) 
where F is volume flow in the unit of m*/h. We can get the volume V by injecting water to the empty 


dummy load, assume the volume flow of water is Fo, the spend time from beginning to inject the water 


to the time when the outlet flow out water is dto. So, 


v=o (32) 
_ Foto 
ôt = F (33) 


In theory, the power can be measured by formula (9) and (33) when the inlet water temperature is variable. 
But in fact, the real-time measurement of the water temperature is difficult. In our system, the maximum 
response time of the temperature transmitter is about 18s. So, this measurement method is not so accurate 
because the temperature measurement delay. Usually, we measure the power by keeping the inlet water 


temperature as a constant. But for verifying this measurement theory is correct, we do some experiments. 


As shown in Fig. 7, the inlet water temperature increases gradually at a constant slope. We can get the 
correct power using formula (9). The definite integral is calculated as the green rectangle shown. The 
temperature difference at time t should be the outlet water temperature at time t minus the inlet water 
temperature at time t-ôt. But if we use formula (10) to calculate the power, the definite integral is 
calculated as the brown rectangle shown, the temperature difference is AT less than the correct value in 


every moment. So, the measurement result using formula (10) will be smaller than the correct result. 


Temperature /°C 
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FIG. 7. The temperature of the outlet water and the inlet water along with the time when the inlet water temperature 
is gradually increasing. The black curve indicates the outlet water temperature and the blue curve indicates the inlet 
water temperature. 


We make the inlet water temperature increase gradually in the gyrotron pulse, and use formula (10) to 
calculate the power. Some test results are shown in Table 2. Especially, shot number 1 is measured when 
the inlet water temperature is gradually increasing in the pulse duration. Obviously, the power test results 


(506kW, 812 kW) are smaller than the real value (546 kW, 849 kW). The temperature difference of the 


dummy load circuit is shown in Fig. 8. The inlet water temperature rose 0.09 °C after 60 s. 


When we use formula (10) to calculate the power, the temperature difference is AT less than the correct 


value in every moment, and, 


ôt =~" x 60 (34) 


~ 0.09 


By formula (10), we can get, 


Fy'GAT‘60 _ 42x4180xATX69 L 2926 x 103- AT = 546 x 103 — 506 X10? = 40x103 (35) 


3.6°T 3.6X1 


where, Fy is the volume flow with the unit of m*/h, change it to the mass flow F with the unit of kg/s, 


_ F:1000 _ Fy 
F = no ae 8/9) (36) 
So, 
AT = 0.0137 °C (37) 
st = x 60 ~ 9.15 (38) 
0.09 


Similarly, for the collector power, we have, 


FV CATO — 391x4180xATX60 — 2724 x 103 - AT = 849 x 10? — 812 x 103 = 37x10? (39) 


3.6°T 3.6X1 


So, 
AT = 0.0136 °C (40) 
Because the inlet water temperature of collector is same as that of dummy load, the inlet water 


temperature rose 0.09 °C after 60 s, 


ôt =Æ x 60 = 9.15 (41) 
0.09 


The analysis above is not very accurate because the temperature measurement is not exactly accurate, 
and the inlet water temperature increased gradually at just a nearly constant slope, but not a real constant 


slope. 


TABLE 2. Power measurement results when the inlet water temperature is gradually increasing (No. 1), and results 
when the inlet water temperature is a constant (No. 2 ~ No. 6). All these data are measured when the pulse duration 
is 1s with the cathode voltage is -44 kV, the anode voltage is 21 kV, the main magnetic current is 50.4 A, the filament 
voltage is 33.4 V, and the filament current is 34 A. 


Average Average 


Dummy Load Probable Collector . Probable 
No. Power / kW Error / % Donny Load Power / kW Collector Error / % 
ý Power / kW Power / kW 2 
1 506 506 812 812 
2 544 847 
3 531 841 
4 560 3.5 546 855 849 3.6 
5 538 849 
6 557 851 
18.2 
18 
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FIG. 8. The temperature of the dummy load outlet water and the inlet water along with the time in shot number 1. 
The blue curve indicates the outlet water temperature and the red curve indicates the inlet water temperature. 


As shown in Fig. 9, when the inlet water temperature is decreasing gradually, the analysis is similar 
with that when the inlet water temperature increased gradually. We can get the correct power use formula 
(9), the definite integral is calculated as the green rectangle shown. But if we use formula (10) to calculate 
the power, the definite integral is calculated as the brown rectangle shown, the temperature difference is 
AT bigger than the correct value in every moment. So, the measurement result using formula (10) will be 


bigger than the correct result. 


Temperature /°C 
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FIG. 9. The temperature of the outlet water and the inlet water along with the time when the inlet water temperature 


is gradually decreasing. The black curve indicates the outlet water temperature and the blue curve indicates the inlet 
water temperature. 


IV. TEST RESULTS 


We have tested our gyrotrons and got the output power using calorimetric method. Up to now, the #1 
Gycom gyrotron has demonstrated 900 kW output power (813 kW power is absorbed by dummy load) 
for 10 s with cathode voltage of -45.5 kV and anode voltage of 24 kV at about 48% electrical efficiency. 
Some related key signals are shown in Fig. 10. And the #1 Gycom gyrotron has demonstrated 644 kW 
output power (585 kW power is absorbed by dummy load) for 754 s[1] with cathode voltage of -44.5 kV 


and anode voltage of 24 kV at about 53% electrical efficiency. 


The #2 CPI gyrotron has got about 499 kW output power (449 kW power is absorbed by aluminum 
tank dummy load and waveguide dummy load) for 79.66 s with cathode voltage of -57 kV and anode 
voltage of 21 kV at about 27.1% electrical efficiency and has got about 721 kW output power (649 kW 


power is absorbed by aluminum tank dummy load and waveguide dummy load) for 517 ms with cathode 
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voltage of -59 kV and anode voltage of 22 kV at about 35.2% electrical efficiency. Some related signals 
are shown in Fig. 11 and Fig. 12. The long-pulse and high-power test is stopped by dummy load leaking 


problem. 
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FIG. 10. The key signals of Gycom gyrotron at 900 kW output power. 
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FIG. 11. The key signals of CPI gyrotron at 499 kW output power. 
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FIG. 12. The key signals of CPI gyrotron at 721 kW output power. 

Except the dummy load and the collector, we measured the power of the other cooling circuits, such 
as Main Window, MOU etc.[15] The test results for #1 Gycom gyroton are shown in Table 3, and the 
probable calorimetric measurement error, water flow, and pressure drop for each cooling circuit are 
shown in Table 4. The power consumption is the sum of the cathode power and the anode power which 
are the product of voltage and current. As we can see, the sum of calorimetric power in all measured 
cooling circuits is approximately equal to the power consumption. It can prove that the calorimetric 


method is self-consistent, and the calorimetric measurement results are credible. 


TABLE 3. The power of all measured cooling circuits for Gycom gyrotron. The power is in the unit of kW. 


No. Pulse Power Dummy Collect MainW Relief Relief Test MOU Output Sum of 
Duratio Consump Load or indow Window Load Bench Power RF Calorim 
n (s) tion Power Power Power Power Power Power Power etric 
Power 
1 99. 961 1005 306 644 0. 62 1.07 4.5 5. 45 23.2 334. 65 984. 84 
99. 981 1114.7 388 652. 4 0. 69 1. 62 6. 56 8. 33 30.32 426.65 1087. 92 


TABLE 4. The probable calorimetric measurement error, water flow, and pressure drop for each cooling circuit. 


Cooling Cuicuit Dummy Collector MainWindo ReliefWindow Relief Load Test Bench MO 
Load wW U 

Probable 

Calorimetric 3.3 3.4 5.2 5.4 2.8 2.8 3.8 

Measurement 

Error (%) 

Water Flow (m?/h) 42 39.6 2:3 2.2 5.2 5.1 3.4 

Pressure Drop 0.3 03 03 0.5 0.3 0.3 0.5 


(MPa) 


V.CONCLUSION 


Calorimetric measurement system has been established to measure the rf power generated by gyrotrons. 
The calorimetric measurements when the water temperature at the inlet is invariable or variable are 
discussed respectively. For getting more accurate results, we usually measure the dummy load power by 
keeping the inlet water temperature as a nearly constant. And the calorimetric measurement uncertainty 
has been discussed. The test results show that the calorimetric method is self-consistent and is credible 


for power measurements. 
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